Multicomponent fluids of hard hyperspheres in odd dimensions 
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Mixtures of hard hyperspheres in odd space dimensionalities are studied with an analytical ap- 
proximation method. This technique is based on the so-called rational function approximation and 
provides a procedure for evaluating equations of state, structure factors, radial distribution func- 
tions, and direct correlations functions of additive mixtures of hard hyperspheres with any number 
of components and in arbitrary odd-dimension space. The method gives the exact solution of the 
Ornstein-Zernike equation coupled with the Percus-Yevick closure, thus extending the solution for 
hard-sphere mixtures [J. L. Lebowitz, Phys. Rev. 133, 895 (1964)] to arbitrary odd dimension. Ex- 
plicit evaluations for binary mixtures in five dimensions are performed. The results are compared 
with computer simulations and a good agreement is found. 
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INTRODUCTION 



Systems made of hard bodies, i.e., impenetrable parti- 
cles interacting solely through hard-core repulsions con- 
stitute useful simple fluid models [l|. Multicomponent 
hard-sphere mixtures serve as important reference sys- 
tems in condensed matter and are relevant for under- 
stating the behavior of complex fluids with additional 
inter-particle interactions, such as those in colloidal sys- 
tems. These simple models capture the main features of 
the packing effects at short distances as they occur also in 
fluids governed by additional, attractive interactions. Be- 
sides, the knowledge of the structural properties of hard 
sphere fluids is a prerequisite for treating attractive in- 
teractions perturbatively and within a density functional 
theory approach. 

The study of d-dimensional hard-sphere fluids may 
prove to be a useful guide in investigating the solu- 
tion of the three-dimensional (3D) problem, apart from 
its own importance at a fundamental level. This ex- 
plains the continued interest on hard-hypersphere sys- 
tems found in the literature along the years [2H74j . How- 
ever, multicomponent fluids at dimension d > 3 have re- 
ceived much less attention and the available information 
is rather sparse. The performed studies reduce to some 
evaluations of equations of state and virial coefficients 
flU l30l . [HI, E3| , phase transition analyses [37| , and com- 
puter simulations [ifjf . To our knowledge, none of the 
widely used mechanical-statistical theories (e.g., integral 
equation theories) have been applied to these systems. 
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The Percus-Yevick (PY) theory [75[ is one of the clas- 
sical approximations of liquid-state theory and, certainly, 
one of the most widely used. In the case of hard particles, 
exact PY solutions have been found for sin gle-c omponent 
fluids in odd dimensions, d = 1 HS, d = 3 j77ll7ci|. d = 5 
S0,d = 7[Il[6l,d = 11 [gg, and also recently for 
even dimensions, d = 2 [79[ and d = 4,6,8 [67}. Never- 
theless, in the case of mixtures, PY solutions were pro- 
vided only for the 3D fluid [§3, EH , apart from the exact 
solution known for the mixture of hard rods [82l |. The 
present work is an attempt to cover the gap on hard par- 
ticle mixtures at dimensions higher than three using the 
so-called rational function approximation (RFA). 

The RFA approach was originally developed for hard- 
sphere fluids j83|. The method was successfully applied 
to other related systems [84[, such as hard-sphere mix- 
tur es I8ll , sticky hard spher es I85rl87| , square- well fluids 
[8cl l89l|. penetrable spheres |9Cll|. and one-c omp onent hy- 
perspheres [fH H3]. As we showed in Ref. J65[, the RFA 
method, in its simplest version, recovers the exact PY so- 
lution for one-component hypcrspherc fluids in any space 
of odd dimension. 

The aim of this paper is to extend the RFA theory 
to additive mixtures of hyperspheres in odd-dimensional 
Euclidean space. It is shown that the RFA method yields 
the exact solution of the Ornstein-Zernike (OZ) equation 
with the PY closure. While the method is generalized to 
any odd dimension, we focus in particular on the solution 
for the five-dimensional (5D) system.. For this system, 
we analyze some of its thermodynamic and structural 
properties in the case of binary mixtures, and compare 
them with available simulation data. 

The remainder of this paper is organized as follows. 
Section |TT] briefly describes some basic quantities of the 
equilibrium theory for multicomponent fluids of hard hy- 
perspheres. In Sec. IIII1 we introduce a Laplace functional 
associated with the fluid structure factors and derive its 
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general properties. In Sec. IIV| we present the extension 
of the RFA method to mixtures. Explicit formulation for 
the 3D and 5D cases are given in Sec. [V] where com- 
pelling arguments about the equivalence to the PY solu- 
tion are offered. Detailed evaluations for binary mixtures 
in d = 5 are presented in Sec. IVI1 Section IVlII is devoted 
to the conclusions. The most technical aspects of the 
paper are relegated to Appendices. 



II. GENERAL BACKGROUND 

It is useful to give here some definitions that will be 
used in the following. Let p be the total number density 
of an A/"-componcnt hypcrsphcrc mixture, let {xi} (with 
i = 1, . . . ,Af) be the set of mole fractions, and let {<x;} 
be the set of diameters. The overall packing fraction is 



i] = J2i=i Vh where 



rji = v d px i a l 



(2.1) 



is the partial packing fraction due to species i. Here Vd 
is the volume of a d-dimensional sphere of unit diameter. 
For d = odd, 



Vd = 



(tt/2) 



(d-l)/2 



d\\ 



(2.2) 



The structure factor Sij{k) for the particle pair 
is given by 



where 



Sij(k) — Xi&ij + pXiXjhij 



h %3 (k) = / dThij(r)e ik 



(2.3) 



(2.4) 



is the Fourier transform of the total pair correlation func- 
tion hij(r), related to the radial distribution function 
(rdf) gij(r) by 

hij(r)=gij(r)-l. (2.5) 
The compressibility factor Z may be written as 

z = -A= = i + £ ^4^(4)' (2 - 6) 



pk B T 



/'</ 



where p is the pressure, ks is the Boltzmann constant, 
T is the temperature, 



Prn = E Xldf 



(2.7) 



In Eq. (|2.6p we have used r/ — VdPPd- 

The isothermal susceptibility \ is given by 



1 (dp 



^-p^^O), (2.9) 
k B T \d P J T{xi} ^ 

where (k) is the Fourier transform of the direct correla- 
tion function Cij(r), which is defined by the OZ equation, 



hij(k) =c lj (k) + p'^xthu(k)ctj(k). 



(2.10) 



i=i 



In matrix form, the OZ relation can be rewritten as (see, 
e.g., Ref. ii|) 



\-c(k)= \ + h(k) 



(2.11) 



where I is the Af x Af unit matrix and c(fc) and h(fc) 
are Af x Af matrices with elements p^XiXjCij(k) and 

Py/XiXjhij^k), respectively. The compressibility equation 
of state (|2.9|) can be written as 

N 

x- 1 = £ v^7[i-c(o)] 

.V 



l + h(0) 



(2.12) 



where in the last step use has been made of Eq. (|2.11[) . 
In particular, in the case of binary mixtures (Af — 2). x 
takes the form 

= [1 + psifen(0)][l + px 2 fe 22 (0)] - p 2 x 1 x 2 h 2 12 (Q) 
1 + p Xl xSiM + MO) - 2h 12 {0)} 

(2.13) 

The zero wavenumber value of hij(k) can be expressed 



a.s 



where 



(2.14) 



dr /iij(r)r" 



(2.15) 



is the mth moment of hij(r). 

III. THE LAPLACE FUNCTIONAL Gij(s) 
A. Definition 



denote the moments of the diameter distribution, and 
<7ij(4) is the contact value of the rdf, cry being the min- 
imum possible distance between particles i and j. For 
additive mixtures considered here one has 



(2. 



In analogy to the case of one-component fluids [65( , we 
introduce the Laplace functional of the rdf in a Euclidean 
space of odd dimension d as 



G ij( s ) = / drrg i3 (r)9 n (sr)e~ 
Jo 



(3.1) 
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which is defined in terms of the reverse Bessel polynomial 
8 n (t) of degree n= (d — 3)/2: 



f=0 



(2n - ■ 



(3.2) 



More details of these polynomials and their properties 
can be found in Ref. [65|. Here, we recall that the Fourier 
transform of the total correlation functions can be ex- 
pressed in terms of Gij(s) as (see Appendix |A"| 



hij(k) = v d 



r.d-2 



s—ik 



where 



-2tt) 



(d-l)/2 



(3.3) 



(3.4) 



The structure factors readily follow from Eq. (|2.3|) . We 
note that the knowledge of Gij(s) allows us to obtain all 
the structural and thermodynamic properties of a multi- 
component hard e?-sphcrc fluid. 



B. Gij(s) at long and short wavenumber 

Being gij(r < cry) = for hard-hypersphere fluids, 
from Eq. (|3.1j) . one obtains at long s [65| 



lim s ^-^e^ s G ij {s) = a^' 2 9ij {ai j ). (3.5) 

This equation determines the contact values of the radial 
distribution function and we will use this to obtain the 
virial equation of the state of the fluid through Eq. (|2.6|) . 

On the other hand, from Eqs. (|3.1[) and (|A6[) . Gij(s) 
may be written in terms of the total correlation function 
as 



Gij(s) 



(d-2)U 



drrhij(r)9 n (sr)e~ 



(3.6) 



As a consistency test, note that Eq. (|2.14j) is reobtained 
from Eqs. ([33)1 and (j33)l . making use of Eq. (f2~2"j) and 
ctn,2n+i = (— l)' i+1 /(2n + 1)!!. Therefore, the expansion 
of Gij(s) in powers of s allows one to identify Hij^-i 
from Eq. Q3.9p . This in turn gives the compressibility 
equation of state via Eqs. (|2.12[) and (|2.14[) . 

In summary, the behaviors of Gij (s) at long and short 
s are directly connected to the thermodynamic variables 
Z and Xj respectively. 



C. Gij(s) at low density 

The lowest terms in the density expansion of the rdf, 
9i j (r)=9%\r)+pgS\r) + 0(p 2 ), (3.10) 



are 



<4 0) (r) = 6(r-^-), 



(3.11) 



9$\r) = 0(r- a tj ) a**W i£ (r), (3-12) 
l 

where Q(x) is Heaviside's step function and f2 a! b(r) is the 
intersection volume of two hypcrsphcres of radii a and b 
whose centers are separated by a distance r < a + b. In 
Laplace space, one has 



G ij (s) = GP(s) + pG^(s) + 0(p 2 ) 



with 



e n+x {(Xi j s)e 



^,(1), 



G f g ) (a)=2J^i(a) I 



(3.13) 



(3.14) 



(3.15) 



The Taylor expansion of e sr yields [6 
(d-2)!! 



G&As 



G l3 {s) = 



POO 

(s) = / &rr6 n {sr)n am {r)e- sr . (3.16) 
a n , m Hij tm +is m , (3.7) In Eq (jgj^j) uge hag been made of Eq (jg^ Tq ob _ 



m=Q 



where the numerical coefficients a n , m are given by 

min(n,m) 



E 

£=0 



(-1) 



m-l 



(m - 1)\ 



(3.8) 



The first n coefficients a n .m with m = 2q + 1 = odd 
(q = 0, . . . , n — 1) vanish [65J. Therefore, 



(d-2) 



„2m 



m=0 

oo 

Hij^m+2s m . (3.9) 



tain Gg](a) we first need the overlap volume Sl a .b(r). 
An explicit expression for the latter quantity, valid for 
arbitrary radii and d — odd, is derived in Appendix [B] 
[see Eqs. (|B10[) and (|B16|) ] . Figure [T] shows graphs of the 
scaled intersection volume between two hypcrspheres of 
radii o\ and a% > o\ as a function of the center separa- 
tion r (in rescaled units) for the first six odd-dimensions 
(d = 1,3,. ..,11), and for several size ratios a^j <J\. 

Inserting Eq. (|B10j) into the definition of G^j(s) we 
get 



(27T) 



(d-l)/2 



dr 



(""if ) 

4n+4 



d-3 



t (sr)e" 



(3.17) 
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1: (Color online) Scaled overlap volume 
^ t T 1 ,CT 2 ( r )/^CT 1 ,(T 1 (0) of two hyperspheres with radii a\ 
and (72 > tJi as a function of the center distance r, for 
Euclidean space with odd-dimension between 1 and 11 (lines 
starting from the straight line), and for several values of 
the size ratio <72/0i. The lower panel corresponds to two 
identical hyperspheres. The limit a-z — > oo with o\ fixed 
represents a hypersphere of radius o\ crossing a flat wall. 



where we have taken into account that cry > \<Ju~ (Jje. \ = 

|<7j — Ojl/2. It can be checked that G^j(s) has the fol- 
lowing structure: 



IV. THE RATIONAL FUNCTION 
APPROXIMATION 

In this section, we describe the RFA method to obtain 
the functional Gij(s), which in turn allows us to obtain 
the structural and thermodynamic properties of the hy- 
persphere mixture. The method presented here is, on 
one hand, an extension to mixtures of the one recently 
applied to one-component systems of hyperspheres [65[ 
and, on the other hand, an extension to hyperspheres of 
the one already proposed for mixtures of hard spheres 
(d = 3) 0. 

The approximation we propose consists of assuming 
the following functional form 



Gij(s) 



[LM-B- 1 (*)]„, 



where L(s) and B(s) are Af x Af matrices given by 



loo 



B(.s) = I 



n+l 

£' 

m=0 



n+l 

m=0 



(4.1) 



(4.2) 



(4.3) 



Here L m , < m < n + 1, arc Af x Af matrices whose 
elements may depend on the fluid density, the particle 
diameters, and the mole fractions of the system, but are 
independent of s. Therefore, L(s) has a polynomial de- 
pendence on s of degree n+l = (d— 1)/2. Besides, TO (s) 
are diagonal matrices with elements given by 



[<M*)]« = "dXivf-" 1 -*'- 



^d-m-sia-is), (4.4) 



with 



(3.18) 

where Qiji(s) is a polynomial of degree 3n+4 = (3d— 1) /2 
which can be further decomposed as 



n+l 



Qijl{s) = S d+1 Q l]t {s) - 9 n +l{aUS) ^2 ^n+l,rn(o- 3 £s) T 



2n+3— m 
9=0 



(3.19) 



where Qij((s) is a polynomial of degree n = (d — 3)/2. 
Note that Qije(s) = Qju(s) but Q ije (s) ^ Qju(s). 

It is worth pointing out that all the equations in this 
section and in Sec. [IT] are exact. 



4>m{x) 



E 

U=0 



a 



(4.5) 



The series expansion of B(s) and its asymptotic long-s 
value are given in Appendix [C] The parameter S in Eq. 
(|4.4|) encompasses two different conditions of normaliza- 
tion for B(s), specifically [see Eqs. (|C8|I and (jCllj) ] 



lim B(s) = I 



(5 = 0), lim B(s) = I (6=1). (4.6) 



These cases constitute two alternative choices of the an- 
alytical representation for Gij(s) which, however, yield 
identical physical results (each one with its particular 
quantities L m ). 

In the functional form (|4.1[) for Gij(s) we are using 
the constraints derived in Sections MI Bl and MI CI In 
particular, as shown in Appendix [Dl Eq. (|4.1[) is consis- 
tent with the exact low-density expansion given by Eqs. 
(I3~T31 . (13~T41 . (1335)1 . (133511 . and (13391 . 

The number n + 2 of terms in the representations of 
L(s) and B(s), Eqs. (|4.2[) and (|4.3[) . is the minimum one 



5 



required to verify the correct behavior of Gij (s) at large 
s. In fact, Eqs. (|33|) and (|4~T|) yield 



(4.7) 



= lip Us) ■ B-\s)l 

= [U+x-B-Hoo)].., 



with B(oo) given by Eq. (|C1 1|) . With the normalization 
choice <5 = 1, the contact values gij{cr^j) are directly re- 
lated to the components of L„ +1 , 

4" +1) ^(4) = (Ln+l)« ( S=1 )- ( 4 - 8 ) 



Now we want to determine the n + 1 coefficients L m 
This is done by requiring consistency with Eq. (|3.9|) 
First, let us rewrite Eq. (|3.7p as 



— G u ( s ) = i + ^(r m ), y . s 



m-2 



(2n + l) 



where we have introduced the matrices V m as 



(4.9) 



(r m )« = "r^:: 1 . (4.10) 



Next, we note that 



(2n+ 1)!! 



where 

Consequently, 

D 2 



(4.11) 



(4.12) 



(2n + l)l! 



e'«'Gf tf («)= J2 ( K ™M"\ (4.13) 



with 



m-2 



K m = A m + ^]r«® A m -e-2, (4.14) 

where the symbol £g> denotes a matrix product "element 
to element": (A <S) B)ij = AijBij. 

Apart from the introduced notation, Eq. (|4.13[) is to- 
tally equivalent to Eq. (|3.7[) . Now, according to the RFA 
form (14.11). 



E L m s m - E K ™ s ™ • E Bfc s " > ( 4 - 15 ) 

m=0 \m=0 / \fc=0 / 

where 



and the matrices are given by Eqs. (|C5[) and (|C6I) 
From a power analysis of (|4.15|) . one obtains 



E ^- fe ' Bfc = U 



(4.17) 



with the convention Le = if £ > n + 1. If we choose 
<5 = 0, the first relation in (|4.17l) (i.e., £ = 0) is trivially 
solved and yields 



(U,)y = (2n+1)!!, (5 = 0). 



(4.18) 



Since, as said below Eq. (|3.8[) . the first n coefficient 
a n m with m = odd < In — 1 vanish, we have 



i~i = r 3 = • • • = r 2 „_! = o. 



(4.19) 



The property (|4.19|) can be used, together with Eq. (|4.17[) 
with £ = even < In + 2, to express the matrices r2,„, 
< m < n, in terms of the matrices \-i by means of the 
recursion relation 

m+l 2k 

r 2 m = Um+2 "EE [ r 2(m~fc) ® A 2fc-^] ' B *, I 4 - 20 ) 

fc=l f=0 

where we must adopt the convention (l~_2)^ = 1. Simi- 
larly, from Eq. (|4.17[) with £ = odd < 2n + 1 one obtains 

m+l 2fe-l 

Um+l = E E t r 2(™-fc) ® A 2 fe_£_l] ■ Bf, (4.21) 
fe=l £=0 

with < m < n. Once the matrices V2m are obtained in 
terms of via Eq. (|4.20[) . Eq. (|4.2ip becomes a closed 
system of n+ 1 = (d— l)/2 matricial algebraic equations 
for the unknowns L2,. . . , L n+ i. Therefore, the prob- 
lem of finding Gij(s) as given by Eq. (|4.ip is reduced to 
solving Eqs. (|4.20[) and (|4.2ip . In the case of hard-sphere 
mixtures (d = 3) it is possible to find an analytical solu- 
tion, as proved in Ref. [8l| . 

The solutions of Eq. (|4.2ip that are physically mean- 
ingful arc those that verify the correct behavior in the 
limit p -> 0. For <5 = 0, we find from Eqs. ([3TT3J) . (j3~14l) . 
and (gUD 



lim (L m ). 



(4.22) 



where the coefficients w n +i m of the reverse Bessel poly- 
nomial n +i(t) are given by Eq. (|3.2p . Alternatively, 
when all particle diameters adopt the same value (say 
Gij = 1 Vi, j), we must recover the pure-fluid solution, 
which only depends on density or packing fraction (but 
not on the mole fractions), i.e., 



(4.23) 



where the functions a m (rj) are defined in Ref. [65|. Ap- 
pendix [E] shows that the one-component solution [(35| is 
actually recovered from the approach (|4.ip . 
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FIG. 2: (Color online) Matrix element of L2 for a binary mix- 
ture at d — 5 as functions of the diameter 02 (crj = 1) for 
X2 = \ and different values of the packing fraction. From top 
to bottom the curves correspond to 77 = 0, 0.01, 0.02, 0.05, 0.1 
and 0.15. The pure-fluid solutions are indicated by symbols. 



Figure [2] illustrates the elements of the matrix L 2 , as 
obtained from the physical solution of Eq. (|4.21[) , for a 
5D binary mixture (d = 5, n = 1) as functions of the 
diameter 02 (with o\ = 1) for xi = j and several density 
values. The dashed lines correspond to the limit of zero 
density, as given by Eq. ()4.22|) . while the symbols show 
the solution for a one-component fluid [see Eq. (14.231) ]. 
with 02(77) given by Eq. (E10) of Ref. (65j|. 

Once the solutions of Eq. (|4.2ip are obtained, they may 
be used in Eq. (|4.7[) to obtain the contact values of the 
pair radial distribution functions and, subsequently, the 
compressibility factor in the so-called virial route using 
the equation of state (|2.6[) . Its expression reads 

(4.24) 

Furthermore, the isothermal susceptibility \ of the 
fluid may be evaluated by means of the thermodynam- 
ical formula (|2.9[) . using the relation resulting from Eq. 
(gHH) and the definition of T m in Eq. (|4~TU|) . 

h(0) = 2(-2^ d - 1 ^ 2 (d - 2)!!r 2n+1 , (4.25) 

with f~2n+i given by 

ra+2 2fc-l 

r 2 n+i = - [ r 2(»-fe+i) ® Aa*-*-i] • Bj, (4.26) 

fe=i £=0 



which proceeds from Eq. (|4.17|) at I = 2n + 3, taking into 
account that \-2 n +3 = 0. The values of the isothermal 
susceptibility can be used then to obtain the compress- 
ibility factor in the so-called compressibility route as 

Z C ( V ) = f dxx-^rjx). (4.27) 
Jo 

Of course, the structure factors and total correlation 
functions can be obtained easily with Eqs. (|2.3[) and p.3[) . 
once the solutions L m have been substituted into Eq. 
(|4.ip . Notice that, as shown in Appendix |D| Eq. (|4.1j) 
is exact to first order in p, so that the first three terms of 
an expansion of the right-hand side of Eq. (|2.3[) in powers 
of p are exact within the RFA method. 

V. EXPLICIT EXPRESSIONS 

In this section, we briefly revise the RFA solution for 
mixtures of hard spheres [d = 3, n = 0) and derive re- 
sults for 5D hyperspheres (d = 5, n = 1). The 3D anal- 
ysis from the present framework is useful as a guide for 
subsequent applications of the RFA method to mixtures 
in higher dimensions. Here we use the version of the RFA 
approach based on the choice <5 = 0. 

For simplicity, henceforth, we introduce the following 
matrix notation 

[Ai} a p = A a) [A } \ aP = A^ 

\A k } aP ee A a 5 a f3, lAijjafi = A a p. (5.1) 

This means that, given a list A = {A\, A2, . . . , Aw}j 
[Aij represents the Af x J\f square matrix made by re- 
peating the list A as columns, so that all the elements 
of a given row are equal. Analogously, [A,] = \A$ 
(where f indicates the transpose of a matrix) is ob- 
tained by repeating the list A as rows, so that all the 
elements of a given column are equal. If the elements 
of A are placed along the main diagonal one gets the 
diagonal matrix \Ak\- The meaning of {A^j is self- 
evident. We will also use the notation [1] to refer to 
a matrix with all the elements equal to 1. Note the 
properties [C fc ] • {Ajj = \C l A j \, \C k \ ■ {A t \ = {C^}, 
I^-J ' ICk] - \C 3 A 3 l {Aij ■ [C fc ] = IMCil Some in- 
version properties involving these matrices are proved in 
Appendix |FJ 

A. Hard-sphere mixtures 

For d = 3 (n — 0), Eq. (|4.ip becomes 

G ^ s ) = e ~^T- [(Lo + Us) • B- X ( S )],. (5.2) 

with 

B(.s) = l+/,[0 o (s)-Lo + < l ) i(s)-Li], (5.3) 



7 



where l_o and Li are the unknowns. From Eq. (|4.18[) we 
have L = [1] and (|43T|) yields 



¥4 



(| + M ). Ll = K-] + ^[l] 



(5.4) 



with /j, m defined by Eq. (pTF)) . Using Eq. ((FT]), it is 
straightforward to obtain 



Li 



M , hi] 



2 2(1 + 27?) 2^ 3 (1 + 277) 



[!]• (5-5) 



This closes the determination of Gy(s). 

In order to evaluate the contact values of gij [see Eq. 

63)], 



°H 9ij(<r?j) = [ L i ' B . 



(5.6) 



we must determine B(oo). From Eq. (|C11|) (with 5 = 0), 
B(oo) = I + [A,] + [Q^J, (5.7) 

with 



_ 3?/ /i 4 77^ 

Al = T - ^ ^ 

1 + 277 /i 3 di 



3 ?7i 



1 + 2?7 cr. 
Using the relation (|F4[) . it results 

B- 1 (oo)= = l + -{apIi]-(H-a)[fiiO-j] 
7 



(5.8) 



(5.9) 



-(l + ^IAJ+ftlA^-]}, (5.10) 



where 



a ee tr([Ai]) = -77 1 + /3 



1 1 ' 



6=Etr(in i ])=/3 



/_*2 



(5.11) 



(5.12) 



tr([A l0i -]) = -2^(1 + 0), (5.13) 



M3 



/3 ee trd^-J) 



7 ee (l + a)(l + /3)-ah 



377 



1 + 27?' 



(I-77) 2 



(5.14) 



(5.15) 



I + 277 

The matrix products on the right-hand side of Eq. ([5 



M-IA,] = al^l [ CTj ].[A 2 ]=aIl], 

• [A^-] = a[a^-], [a,] • [A^-] = a fa], 
[1] • [A,] = 41], [1] • [Ai^-J = afa], (5.16) 



plus similar results in the case of [fij]. Taking into ac- 
count the preceding results, one obtains from Eq. (|5.6[) 



1 377 M 2 o"»o-j 

1 - 77 2(1 - 7/) 2 /J 3 (Ty 



(5.17) 



Its application in Eq. (|2.6j) gives the pressure equation in 
the virial route obtained by Lebowitz [80| . 

On the other hand, thanks to Eq. (|5.5|) . it is possible 
to solve Eqs. (|4.25j) and (|4.26[) analytically. The result 



is 



M°) »7 /9t7 2 (/j 4 //J3) 3 6t7/j 4 /j 5 //x§ /j 6 



tt 8 V (1 + 277)2 I + 277 ' Ms J ^ 

(1 + 2t7)/j 5 //J3 - 3t7(/j 4 //J3) 2 . , 1, 



(1 + 277)2 



II O? + tr?l 



8 L 



^Wm3)(I4I + k^j-I) - jgfa + ogjj 

2(1 + 2t?) 



(5.18) 



The use of (f5~T8"l) into (|2~T3"1) yields the PY solution for 
the isothermal susceptibility Xi which agrees with that 
given by Ashcroft and Langreth (9l| for a binary mixture 
of hard spheres. 



B. 5D mixtures 



For a fluid of hyperspheres in d = 5 {n = 1), the func- 
tional Gij(s) takes the form, 



with 



[(Lo + Ua + lVj-B- 1 ^, (5.19) 



B(a) = \ + p[0 o (s) • L + «t>i(«) • Lj + * 2 ( s ) • L 2 ] . (5.20) 



According to Eq. (|4~Tg)l , L = 3[1]. Besides, (|4~2"Tj) at 
?7i = yields |92| 

Li = SK-J+y"^ ( 2 VM + 3771a,] - 10fo/^] • L 2 

(5.21) 

This expresses l_i in terms of L 2 . From Eq. (|4.20[) . with 
m = 0, we find 



.'9/J7//J5 4t7(/j 6 //j 5 ) 2 \ 277(/j 6 /// 5 )[cT i jI 

I -^rr- + 5 - ) w rr^ — 

(5.22) 



14 I-677 
"2(1-6??) 2 1 yJ 



Qa- U 



where 



Qa s I 3 | j + 20 ^[^M-I + (5 . 23) 
3 /j 5 (1-6t7) 1 - 677 



Finally, using Eq. (|4.21[) with m = 1, and the relations 
(|5.21j) and (|5 . 22[) , we obtain the following quadratic equa- 
tion for L 2 : 

= Q + Qi • L 2 + Q 2 • L 2 • (P + Pi • L 2 ) , (5.24) 



8 



where 



2M + # + ^M, (5.25) 



Pi = - 



60 



1 - 677 



l-6rj ^5(1 — 677) ' 



(5.26) 



Qo 



■v- 



20r] 2 (fi 6 /fi 5 ) 3 _ 4r)n 6 n 7 /nl 

1 — 6?7 8/J5 
(1 - 677)^7/^5 i, ti 



(1 - 6^)2 

1077(/j 6 //a 5 ) 2 - 



(1 — 6T?) 2 



3r? 



Ms(l - 677) L 
(1 + 24t?)m6 



4/z 5 (l - 6tj): 



J 2(1 



677) 



l<Ji<Jj(Ji 



(5.27) 



Qi 



1077/i 6 //z 5 



foil 



1 - 6?/ r,JJI ' 6 2(1 - 677) 

10077(/V M5 ) 2 + 10(1 - 677) M7 //i 5 
^ (1 - 677)2 

25t7M 6 /M5 it_ , ti 5|gf^/Pj] 



(1 - 677) 2 



2(1 - 677) 



2 

(5.28) 



Due to the fact that Eq. ()5.24j) is quadratic in l_2 the 
evaluation of the structure functions is now considerably 
more complex than in the three-dimensional case. In gen- 
eral, solutions of (|5.24[) must be worked out numerically. 

Binary systems can be completely specified by the total 
packing fraction 77, the concentration of one component 
(say X2), and the diameter ratio 02 /ci- For arbitrary 
and finite values of these parameters, Eq. f|5.24[) yields 
four solutions L2 that can be obtained analytically, one 
of which being the physical root that verifies the conver- 
gence conditions P~2"2l and P~2"3"|) . 



C. Direct correlation functions 

The knowledge of the Laplace functions Gij(s) allows 
one to obtain the direct correlation functions via Eqs. 
(12. lip and (|3.3[) . Although we have not attempted a for- 
mal proof, we have checked that in the binary case the 



structural form of (kj(k) is 

%( fc ) = ^[P^me^+V^i-i^e-^ 
+ Qy(ifc)e i(CT '-^ )fc/2 + Qiji-iky^ 



,)fe/2 



(5.29) 



where Vij(s) and Qij(s) are polynomials of degrees 
(3d — l)/2 and d — 1, respectively. Moreover, the quan- 
tity enclosed by square brackets is of order k 2d , so that 
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FIG. 3: (Color online) Compressibility factor as a function 
of the packing fraction r\ for a binary mixture at d = 5 with 
02/cri = j and mole fractions X2 
on the plot) 

the compressibility routes (solid lines on the right and left of 
the shaded area, respectively) are compared with predictions 
of an analytical equation proposed in Ref. 34[] (short dashed 
lines) and simulation results from Ref. [40] (symbols). The 
long dashed lines represent the interpolation Z = %Z V + \Z C . 



,(0) is finite, as required by Eq. (|2.12[) and enforced a diameter ratio <72/0i 



through Eq. (|3.9[) . The direct correlation functions in 
real space, Cy (r) are obtained from Eq. (|5.29[) by means 
of Eq. (|A2|) . The important point is that, upon applica- 
tion of the residue theorem, the structure given by Eq. 
(|5.29|) implies that Cjj(r) = for r > Oij [65[. Since 
the RFA also complies with the physical requirement 
9ij( r ) = for r < dij, we recover the two conditions 
defining precisely the PY closure to solve the OZ equa- 
tion. Therefore, we find that the RFA developed in this 
paper yields the PY solution for hard-hypcrsphere mix- 
tures of odd dimensions. This is a remarkable result since 
both approaches (RFA and PY) are, in principle, rather 
independent. 



VI. RESULTS FOR 5D BINARY MIXTURES 

Here we consider a mixture of two types (i = 1, 2) of 5D 
spheres with arbitrary diameters 0^ and concentrations 
Xi- For simplicity we fix a\ = 1. 

Figure [3] compares the compressibility factors calcu- 
lated from the virial (Z v ) and compressibility (Z c ) routes 
[cf. Eqs. (|4.24p and (|4.27|) . respectively] with those from 
computer simulations (molecular dynamics) |4(j and a 
semiempirical equation of state [3~H ]. for a mixture with 

i and mole fractions 22 = j, 
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FIG. 4: (Color online) Compressibility factor obtained from 
the RFA-PY method for a binary mixture at d = 5, as a func- 
tion of the particle diameter a-2 (<Ji = 1) for r\ = 0.15 and 
mole fractions X2 indicated on the plot (solid lines). Predic- 
tions of the analytical equation from Ref. [34| are shown with 
dashed lines. 



i and |. We observe that Z c and Z v bound both the 
simulation data and the proposal of Ref. [H|, with Z c 
slightly above and Z v slightly below. The agreement is 
very good at low densities (packing fraction lower than 
0.1) and reasonably good over the whole density range 
(0 < 77 < 0.19) of fluid phase predicted for the one- 
component system. In general, the compressibility route 
gives better results than the virial one. A similar relation 
Z v < Z < Z Cl with Z being the actual compressibility 
factor, has been observed for the PY solution of the pure 
5D system [65| . Figure [3] also includes the interpolation 
formula Z = \Z V + \Z C dH. This Carnahan-Star ling- 
like recipe presents a very good agreement with simu- 
lation data and is practically indistinguishable from the 
semi-empirical equation of state proposed in Ref. [34| . 
except for X2 = \ and r\ > 0.15. Unfortunately, no sim- 
ulation data are available in those cases. 

Next, it is instructive to examine the size ratio depen- 
dence of the compressibility factor. Figure 2] compares 
values of Z v obtained from the RFA-PY solution via the 
virial route with those calculated from the analytical ex- 
pression given in Ref. |IH . These evaluations correspond 
to a packing fraction r) = 0.15 and a set of molar fractions 
ranging from xi = ^ to X2 = j- As the concentration X2 
decreases, we find, as expected, that the compressibility 
factor comes to have a weaker dependence on 02 ■ The 
general trends of Z v versus 172 for different mole frac- 
tions X2 qualitatively agree with the predictions of the 
equation of state proposed in Ref. [34j (there are no sim- 
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FIG. 5: (Color online) Number-concentration structure fac- 
tors Snn, Snc and S cc , as obtained from the RFA-PY method 
for a binary mixture in d = 5 dimensions with diameter ra- 
tio 0-2/eri = j, packing fraction r\ = 0.15, and mole fractions 
X2 = \, j, |, yg, and ^ (solid lines). As X2 is reduced, 
the curves converge to the solution of the pure fluid (dashed 
lines). 



ulation data available for this analysis). The differences 
between both calculations (at most 9%) are similar to 
those observed before. 

Now we analyze the structure and correlation func- 
tions. Instead of the conventional structure factors Sij, 
we consider some combinations of them which may be 
easily associated with fluctuations of the thermodynamic 
variables (Ull!: 



S NN (k) = S n (k) + S 2 2(k) + 2S 12 {k), 



3.1a) 



SNc(k) = x 2 Sn(k) - xiS 2 2(k) + (x 2 - xi)Si 2 (k), (6.1b) 



S cc (k) = x 2 2 S n (k) + x\S 2 2{k) - 2x!X 2 5i 2 (fc). (6.1c) 

In the limit of small wavenumber (k — > 0), Snn an d 
S cc become the mean square fluctuations in the particle 
number and concentration, respectively, whereas Sn c is 
the correlation between these two fluctuations. 

We plot in Fig. [5] the number-concentration structure 
factors for a packing fraction r/ = 0.15, a diameter ra- 
tio Oija\ = j, and decreasing values of the particle-2 
concentration, X2 = 5, \, \, jg, and As expected, 
S cc and Snc reduce smoothly to zero with X2 for all 
wavenumber, whereas the number-number structure fac- 
tor Snn converges to the static structure factor of the 
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FIG. 6: (Color online) Direct correlation functions obtained 
from the RFA-PY method for an equimolar (x\ = X2 = \) 
binary mixture in d = 5 dimensions at 77 = 0.15 and various 
diameter ratios, 02/cri = i (dotted line), | (dashed line), 
I (long dashed line), y§ (dashed-dotted line), and || (long 
dashed-dotted line). The pure- fluid solution is displayed for 
comparison (solid line). 



one-component fluid in this limit (dashed line). Notice 
that Snn and Sec are positive for all k by definition. The 
main modifications of the structure factor Snn with in- 
creasing X2 are the variation of height of the main peak, 
with its position very weakly altered, and an increase of 
the values at short k. 

Evaluations of the direct correlation functions defined 
(in the Fourier space) by the OZ equation (|2.10[) are 
shown in Fig. O The results correspond to an equimolar 
mixture for different diameter ratios from 02/ '&i = h to 
nearly equivalent hypersphere sizes, <Jil<j\ = f^. In the 
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limit <t 2 1 o-y — > 1 the three correlation functions Cy (r) be- 
come identical and match the pure- fluid values. Figure [6] 
clearly shows that the direct correlation functions Cij(r) 
vanish for r > <jy , thus confirming that the RFA method 
yields solutions of the PY closure to the OZ equation, as 
discussed in Sec. IV CI 

Figures [7] and [5] show the radial distribution functions 
gij(r) of a highly asymmetric binary 5D-sphere mixture 
with CT2 /(7i = i and x 2 — j, at the reduced densities 
po\ = 1.4 and po\ = 2.0, respectively. The selected 
densities lie in the proximity of the phase transition for 
the pure one-component fluid (77/ = 0.19) predicted by 
Michels and Trappeniers [1]. We see that the RFA-PY 
solution at 77 rj 0.17 gives reasonable good values for 
all pair functions throughout the r- region (see Fig. [7]). 
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FIG. 7: (Color online) Radial distribution functions for a bi- 
nary mixture in d = 5 dimensions with parameters X2 = \, 
o 2 /(Ji = \ and density pa\ = 1.4 (77 = 0.172774). Solid 
lines: results from the RFA-PY method. Symbols: computer 
simulations 14011. 



At a higher packing fraction, 77 m 0.25, we observe in 
Fig. [5] that similar quantitative trends in (722 and 1712 
are followed by both simulations and RFA-PY solutions. 
However, there are discrepancies particularly severe in 
the values of gn corresponding to pairs of big particles. 
The strong oscillations observed in the simulation data 
of (?n can be considered as a signature of a solid phase at 
this density |4(j, which is not captured by the RFA-PY 
solution. 



VII. CONCLUSIONS 

In this paper we have extended the RFA method to 
multicomponent systems of hard-hypcrsphcres in odd di- 
mensions. The main features of this approach are based 
on simple physical considerations on Laplace function- 
als Gij(s) of the radial distribution functions, which are 
closely related to the structure functions of the fluid. The 
basic physical requirements upon which the method is 
based are: (i) The radial distribution functions take fi- 
nite values at contact and vanish inside the core, which 



implies that lim s 



s (5-d)/2 e <r t] s Gij ( s ^ = finitc . 



the isothermal compressibility is finite, which implies 
lim s _».o [Gij (s) — {d — 2)!!/s 2 ] = finite; and (iii) the first 
three terms in the series expansion in powers of density 
of the structure factors are exact. Condition (iii) involves 
the evaluation of the overlap volume of two arbitrary en- 
spheres. For space of odd dimensionality, we have found 
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FIG. 8: (Color online) As in Fig. [Jj but for pa\ = 2.0 (rj = 
0.246820). 



the exact, analytical and closed- form expression for the 
overlap volume of two hypersphcres of arbitrary sizes as a 
function of the center-distance [cf. Eqs. (|B10|) and ()B16[) ] . 
We have been able to perform this evaluation thanks to 
the use of reverse Bessel polynomials and Fourier analy- 
sis. 

The primary result of the present work has been to pro- 
vide a theoretical method for the evaluation of thermody- 
namic and structural quantities of multicomponent mix- 
tures of hard hyperspheres at odd dimensions. We have 
shown that this approach gives the exact solution of the 
OZ equation with the PY closure. From that perspective, 
our work extends, on one hand, the Lebowitz solution 
for 3D-sphere mixtures [8(j to higher dimensions and, 
on the other hand, the PY solution for one-component 
hyperspheres [!, 0, HE [fill to the case of mixtures. 

Although the theory here developed is general and 
applies to any odd-dimensional, hard-hypersphere fluid 
with an arbitrary number of components, the results in 
this paper concentrate on 5D binary mixtures. We have 
checked that the virial (Z v ) and compressibility (Z c ) 
routes to the compressibility factor under- and overes- 
timate, respectively, the simulation data, the interpola- 
tion approximation Z = \Z V + %Z C providing excellent 
results. In addition, reasonable agreement between RFA- 
PY results and simulation data for the radial distribution 
functions are found at densities lower than that of the 
fluid-solid phase transition predicted for a pure fluid at 
this dimension. 

The implementation of the RFA method developed 
here involves as many unknowns as the minimum num- 
ber required by the physical conditions and that is why 



it coincides with the PY solution. However, as done for 
3D mixtures (8i| and for d-dimensional one-component 
systems [65[, one can go beyond the PY level by adding 
an extra matricial term L„ + 2 in Eqs. (|4.2p and (|4.3[) . and 
replacing I by (1 + us)\ in Eq. (|4.3|) . The elements of 
L, i+ 2 and the parameter u are free and can be fixed, 
for instance, by imposing given expressions for the con- 
tact values gij(<r^~j) and the thermodynamically consis- 
tent isothermal susceptibility \- 

We expect that the results presented in this paper can 
stimulate simulation studies on multicomponent systems 
of hard hyperspheres. We also plan to undertake the 
investigation of possible fluid-fluid demixing transitions 
predicted by the equations of state obtained here within 
the RFA-PY approach. 

Acknowledgments 

We would like to thank M. Gonzalcz-Mclchor, who 
kindly sent us her simulation data which were used in 
Figs. [7J and [5] The work of R.D.R. has been supported 
by the Consejo Nacional de Investigaciones Cientificas 
y Tecnicas (CONICET, Argentina) through Grant PIP 
112-200801-01474. A.S. acknowledges support from 
the Ministerio de Ciencia e Innovation (Spain) through 
Grant No. FIS2010-16587, partially financed by FEDER 
(Fondo Europeo de Desarrollo Regional) funds. 

Appendix A: Fourier Transform 

In odd-dimension spaces, the Fourier transform of an 
arbitrary radial function £(r) and the inverse operation 
can be evaluated in terms of the reverse Bessel polyno- 
mial 6 n (t) of degree n = (d — 3)/2 as [§5| 

m = ( > kd _ 2 i J drr£(r)6 n (ikr)e- lkr , (Al) 

e(r) = 1 ' d _ 2 i/ dkkak)0n(ikr)e- lkr . (A2) 

Alternatively, 

£( fc ) = 1 2Im{^ w K(r)](-ifc)}, (A3) 

= 2Im{J n [£(fc)](-ir)}, (A4) 

where T n is a Laplace functional defined by 

^K(r)](«)H / dxxax)0 n (sx)e- sx . (AS) 
Jo 

From an integral relation on 8 n (t) one arrives to [6o| 

Fn[M(s) = — ^ — = ^ , (A6) 
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T n [Q(r - a)](s) 



[as)e 



(A7) 



Appendix B: Co- volume of two hyperspheres 

The overlap volume of two hyperspheres of radii a and 
b. whose centers are a distance r = |r|, can be evaluated 



n„ |t (r) = J dr'e(a - \r'\)Q{b - \r' - r|). (Bl) 

Henceforth, without loss of generality, we assume a < b. 
The right-hand side of Eq. (|B1[) is a convolution. Thus, 
in Fourier space Eq. (|B1[) reads 



(B2) 



(B3) 



fio,*(fc) = ®a(k)e b (k), 

where we have called 

o (jfc) = / dr<d(a-r)e- ik -' 



Now using Eqs. ([A"3j) , (|A6|) , and (|AT|) . and taking into 
account that 0(a — r) = 1 — 0(r — a), we get 



e a (fc) 



( 27r )(rf-l)/2 



21m 



9 n+1 (as)e- as 



(2n+ 1)!! 



s— — ik 



(2tt 



,(d-l)/2 



fc <i 



-i [#„+i(-ifc a )e ifca 



(ifca)e 



(B4) 



Therefore, 



(27T 



fc 2d 



9 n+1 {ika)9 n+1 {-ikb)e^ b - a) 

+9 n+1 (-ika)0 n+1 (ikb)e- ik ^ 
-6 n+1 {-ika)e n+1 {-ikby k{b+a) 



1 (ika)9 n+1 (ikb)e- ik( - b +^ 



(B5) 



We can obtain f2 a .&(r) using the inverse Fourier trans- 
form (|A2|) . Integration on the complex plane and the 
application of the residue theorem yield 



(2 7r )( d - 1 )/ 2 



4;f 4 (r)e(6 + ft -r) 



i?t% 6) (r)e(6-a-r) 



with 



R 



(±a,b) 
4n+4 ' 



r) 



Res 

t=o 



6 n+ i (Tta)6 n+1 (~tb)d n (tr)e 



t(b±a-r) 



fAn+5 



J —^[Ot^0 n+1 ( T ta)0 n+1 (-tb) 



x6 n (tr)e 



t(b±a—r) 



t=0 



(B7) 



Equation ()B7j) implies that J?4„+4 (r) are polynomials 
of degree 4n + 4 = 2c? — 2. If r < 6 — a the smaller 
hypcrsphere will be fully contained within the bigger one, 
so that f2 Q) b(r) will be equivalent to the volume of the 
former, i.e., 



Sl a , b (r) = (2a) d v d 



2a d (27r)( d - 1 >/ 2 
d\\ 



, r<b-a. (B8) 



Thus, according to (|B6|) . the difference between the poly- 

ing cxprcssioi 

2 a 2n+3 r 2n+l 



nomials i? 4 ^ + 4 takes the following expression (for any r) 



^4n'+4 ( r ) ^4n+4 ^ ( r ) 



Therefore, the hnal result is 



(2n + 3)!! 



(B9) 



r ^-^ g d 

d\\ ' 



r < b — a, 



(27r) (d-i)/ 2 ^M > 6 _ a < r < 6 + a) 



0, r>6 + a. 

(BIO) 

Let us now obtain an explicit expression for 
the polynomial i? 4 ^' +4 (r). According to Eq. (|B7|) . 
it is the coefficient of t An+A in the expansion of 
9 n+1 (-ta)8 n+1 (-tb)9 n (tr)e t< ~ b+a - r *> in powers of t. First, 
note that 

2n+2 

e n+1 (-ta)e n+1 (-tb) = cjf **> (Bn) 



e=o 



where 

la,b) 

Next, 



min(^,n-(-l) 
£ 1 =niax(0,^-n-l) 



(B12) 



3)1+2 



9 n+1 (-ta)e n+l (-tb)6 n (tr) = ]T d { n a ' b) (r)t e , (B13) 



£=0 



with 



m'm(£.n) 

E < 

£i=max (0,£-2n-2) 



(B14) 



Finally, the coefficient of t l in the t-cxpansion of 
(B6) 0„ + i(-ia)0„ + i(-t6)0„(ir) e t ( fc + a -'-) is 



E 

£i=max(0,«-3n-2) 



Setting I = 4n + 4 and inserting Eq. (|B14[) . one obtains 

, + ft _ r )3»+2- 

(4n + 4-£)! 



«4(0 = (a + b-r)^j:^±^l 



f=o 
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E UnM^f-t/ 1 - (B16) 

«i=max(0,^-2n-2) 

Thus, for example, at d = 3, 5 (respectively n = 0, 1) one 
finds 



R { t h \r) 



4 a ' b \r) 



(a + b- r) 2 
24 



[r 2 + 2(a + 6)r-3(6-a) 2 ] , 
(B17) 



(a + b-rf r 5 



1920 



[3r 5 + 9(a + 6)r 4 



-2(a 2 - 18ab + b 2 )r 3 - 2(b - a) 2 {a + b)r 2 
+15(6 -a) 4 r + 5(6 - a) 4 (a + &)] . (B18) 

It can be checked that R^%(r) admits the following 



structure: 



Rt%(r) = {a + b-rT+ 2 \r 2 ^P^{r) 



+(b-a) 2 Kt b \r) 



(B19) 



where P%+i{r) and K^n ( r ) are polynomials of degree 
n + 1 and 2n, respectively. Clearly, in the case of iden- 
tical hyperspheres (a = b) R^l^{r) adopts the known 
expression 

Rt%{r) = (2a - r T+ 2 r 2 <^ P%< »(r), (B20) 
with P}W(r) given by Eqs. (B7) and (B8) of Ref. (H. 



with 



Bo = 



I (5 = 0), 

1 + P Elto C d-rn ■ Lm (5=1), 



(C5) 



ra+1 



B* = p E c «W-m • Lm (£>!)■ (C6) 



m=0 

According to Eqs. (|C2j} and ([C4| 
lim0 m (s) = J ^ 

8-K) C d -rr 



(8 = 0), 

(5 = i), 



lim B(s) 



(5 = 0), 



(C7) 



(C8) 



As for the behaviors in the limit s — > oo, it is straight- 
forward from Eqs. P~3" ]) -([3~5 ]l to get 



lim (j> m {x) 



(-iy 



lim * m (s) 



-Q-m (5 = 0), 

(5=1), 



(C9) 



(C10) 



lim B(s) = { ' P Q " m ' Lm ^ ~ °}' (Cll) 

s^oo I I (() — 1). 



Appendix C: Properties of B(s) 

The series expansion of the function <j)(x) defined by 
Eq. (Oil is 



»(*) = -£ 



m+£ 



(m + i)\ 



Therefore, Eq. (jO)) yields 

oo 

<S>m(s)= E C d+e-mS e (5 = 0,1), 



=1-5 



where the diagonal matrices C m are 

-i)/2 T . 



(Cm)« = "(-27T) 



(CI) 



(C2) 



(C3) 



Thus, the series expansion of the matrix B(s) defined by 
Eq. g3J) is 



Appendix D: The Rational Function Approximation 
in the low-density limit 

Let us consider the low-density expansions of L m and 
B(s), 



L^+pLW+fV), 



where 



n+l 



'(*) = E 



m=0 



L(o) 



Insertion into Eq. (|4~Tj) yields Eq. (|3~T3| with 



_n-.., n+l 



E [e 



m=0 



(Dl) 



B(s) = l + pB«(s) + 0(p 2 ), (D2) 



(D3) 



(D4) 



B(s) = ^B^ 

£=0 



(C4) 



— CTyS "+ 1 



E[l«-l#-B«(,)l a-. (D5) 



m=0 
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Comparison between Eqs. (|3.14[) and Eq. (|D4j) implies Using the mathematical property (|F1[) and taking into 
that account that tr([^]) = rj, 



L (o) 



>j 



-n +1 , m <- (D6) B - 1(g) = ,_ A pS^f^L 

1 + Ad?7Lm=o 9d- m {s)a. 



Let us consider now the matrices Lm'' ■ It is convenient 

to express their elements as From E q s - dE2} and (Eg) it readily follows: 

L«l = E^W (° 7 ) L( ' ),B " 1(,) = u>^! m M - (E5) 

Using Eq s gH ) and (@~5]), one can recast Eq. jD5| into Finally, Eq. g3|) becomes 

the form (|3 . 1 5[) with _ s „ m 

G i: j{s) = — ^rn^nS 

m=o This coincides with the PY solution of the one 

n+i component fluid in d = In + 3 dimensions, as derived 



— O n+ i(a ie s) E w„+i, m {cr 3 is) m from the RFA in Ref. [if 



—Vd- 

S" 1 - 

m=0 

d-m-5 , -.g 

x S ' - — — — - (D8) Appendix F: Inversion of matrices 

q=0 9 ' 

Let us consider square matrices defined following the 

Comparison with Eqs. (038} and (|3J9|) allows us to iden- notation of Eq . ^ Therij aas proved beloW) one has 
tify 

n+l 



E L ™!^ s ™ = ^sQ ijt {s) (5 = 0), (D9) 



m=0 



(Fl) 



r(i) m 7T ^ /i / x w »+i,to + lCiBj]) = I - \ J \ F2a) 

E L «4^ S = v d sQ ijt {s) - v d e n+1 {a a s) 2^ 7J—\\ l + P 



(d-m)\ 

rn—0 rn—0 

xa%(-cr i ) d - m ((5 = 1). (D10) 



(I + {BtCj})^ 1 = I - I^i^I (F2b) 
Therefore, the RFA proposal (|4.1[) (|4.5[) is consistent 1 + p 

with the exact result to first order in density. 

(I + [A,] + [C^-ir 1 = \ + -{a{B 3 ]-(l + a){C l B ] ] 
Appendix hj\ Pure-fluid limit 7 

-(l + fll^l + fclQA,-]}, 

We consider here cr^- = <j = IV i,j, within the case 6 = (F3) 
0. Then p4> m (s) = ^d(l>d-m(s)lrikl and one can expect 
that L rn = a r „[l], with a m being a certain constant. Here 

we use the matrix notation defined in Eq. dSTTJ) and have (| + J^.] + [B^-]) -1 = I + - {a[Sj - (1 + a)^^] 
defined 7 

^^-w^, (E1) -a+fflM+^D^ 

Thus, where 

L(s) = E «ms m [l] (E2) a ee tr([^I) = tr([^]), a = tr([a^-]) = trd^Q]), 

(F5) 



and 

B(s) = I + A d E & 

— m (5') ft?n 

foil- (E3) 



6 ee trfl^J) = tr([B 4 ]), /3 ee trQC^]) = trQi^-J), 

n+l (F6) 



m=0 



7= (l + a)(l + /3) -a&. (F7) 
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Equations (|F1[) and (|F2[) follow immediately from the 
properties [A*] • [A,-] = a [A,-], [A] ■ [A] = a[A], 
\BiCj} ■ \BiCil = PiBiCjl and {C.B,} ■ [QB,] = 

The proof of Eqs. (|F3|) and (|F4[) is somewhat more 
involved. We first multiplicate (I + [A,])" 1 by I + [A,-] + 
[dBj] [seeEq. ([FT])]. 



i + Kb| (F8) 



where in the last step we have used that [A,-] • [Cj-Bj] = 
and have introduced 



where we have used (1 + /3')(1 + a) =7. Finally, from 
Eq. (|F9P it is straightforward to see that the right-hand 
side of Eqs. (JK3J) and (|F12|) coincide. 

An analogous method can be used to prove that the 
right-hand side of Eq. (|F3|) is also the right inverse of 
I + [A,] + ICiBjj. Multiplying the latter matrix by (I + 
[aJ)" 1 [see Eq. ([FT])], one gets 



(I + \Aj\ + [dBj]) ■ (I + [A,-])- 1 = I + IdB'X (F13) 



with 



1 + a 



Thus, the left inverse of I + [AJ + {C l B j } is 
(\ + {C[B J \)- 1 -{\ + [A^- 1 -^ {C ' iBj] 
where 



1 + /3' 



p' = tx{\C' i B j \) = p- 



1 + a 



(F9) 



JM 
1 + 

(F10) 



(Fll) 



and use has been made of Eqs. (|F1[) and {F2a|. Taking 
the p roperty {C-Bjj ■ [Aj] = 6[C-Aj] into account, Eq. 
(|F10|) becomes 



l + a 1 + P> 
'-[ClAj], (F12) 



B, = B, Aj. 

3 3 l + a 3 



Thus, the right inverse of I + [A,] + \C l B 3 \ is 



(F14) 



{\ + \A j })- l -{\ + lC i B' j \)- 



VP 
(F15) 



Here we have used /?' = tr([ C t ^ ] and {Ajj 
afBjl. Replacement of Ec;. (|F14|) makes the right-hand 
side of Eq. (|F15[) coincide with the right-hand side of Eq. 
(|F3]) . This completes the proof of Eq. f[F3]l. Equation 
(|F4|) is just the transpose of Eq. (|F3j) . Note finally that 
Eq. jF3]) reduces to Eqs. ((FTJ and (|F2aj) by particulariz- 
ing to (7j = and Aj = 0, respectively. 
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